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Abstract 

Ordered protein layers are the subject of active biomedical research for their 
usually interesting physicochemical properties, e.g. permeability, stiffness and 
pours structure. In presented work, we focused on layers build of fibrinogen 
molecules characterised by strong shape anisotropy. Using Random Sequential 
Adsorption (RSA) method, we simulated adsorption process in which the ori- 
entation of adsorbate was described by a non-uniform probability distribution. 
Thus obtained covering layers had different level of global orientational order- 
ing. This allowed us to find dependence between main properties of layers, such 
as maximal random coverage ratio, and order parameter. For better description 
and deeper understanding of obtained structures, the autocorrelation function 
as well as distribution of uncovered space were determined. Additionally, we cal- 
culated the Available Surface Function (ASF), which is essential for determining 
adsorption kinetics. 
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1. Introduction 

Structure of monolayers formed during adsorption process is strongly related 
to properties of underlying substrate, environmental conditions and adsorbate 
molecules themselves. Earlier works on structured protein coverages indicate 
a high potential of such layers for material, medical and food sciences as well 
as pharmaceuticals and cosmetics industries. In the present study we assumed 
that strongly anisotropic fibrinogen molecules tend to align along specific axis. 
Experiments performed earlier show that ordering can have considerable impact 
on covering layers, e.g. [J, Q. 0m aim is to find out how the ordering affect 
maximal random coverage ratio, coverage structure, and adsorption kinetics. 
To achieve that, we performed extensive numerical simulations on fibrinogen 
adsorption using RSA algorithm. 
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2. Model 



Fibrinogen molecule is modeled as a component of three spheres of 6.7, 5.3 
and 6.7nm in diameter, connected by two chains of ten small l.5nm spheres; see 
FigUl Such a model was used earlier in 0, Q] and it turned out to be the most 



Figure 1: Approximation of fibrinogen molecule shape used in simulations. The 
side balls have diameter of Q.lnm whereas the middle one is a little smaller 
5.3nm. Each of the small spheres between them is of 1.5nm diameter. 

effective in reproducing maximal coverages and adsorption kinetics obtained in 
experiments. 

Fibrinogen molecules are placed on a at homogeneous collector surface ac- 
cording to Random Sequential Adsorption (RSA) algorithm Q, which itera- 
tively repeats the following steps: 

1. a virtual fibrinogen molecule is created with its centre position on a collector 

is chosen randomly according to the uniform probability distribution; 

2. the molecule orientation is chosen randomly according to normal probability 

distribution having specific expected value (p and variance a 2 ; 

3. an overlapping test is performed for previously adsorbed nearest neighbours 

of the virtual molecule. The test checks if surface-to-surface distance be- 
tween each of the spheres is greater than zero; 

4. if there is no overlap the virtual molecule is irreversibly adsorbed and added 

to an existing covering layer. Its position does not change during further 
calculations; 

5. if there is an overlap the virtual fibrinogen molecule is removed and aban- 

doned. 

The number of RSA iterations N is typically expressed in dimensionless time 
units: 

t = N^, (1) 

where Sf = 127.918nm 2 is an area covered by single fibrinogen molecule and 
Sc is a collector size. In case of our simulations algorithm was stopped after 
t = 10 4 . Example coverages obtained this way are presented in FigJ5] To 
check the importance of boundaries shape, simulations were performed using 
both round and square collectors. Obtained results do not show any significant 
dependence on collector shape. For further analysis, coverages generated using 
round collector of 500nm diameter were chosen. 

As mentioned above, the orientational order inside a layer is controlled by 
a parameter a 2 bound up with the width of Gaussian probability density func- 
tion used during simulations. However in possible experiments, ordering can 
be evoked by many different factors, e.g. collector structure, electrostatic in- 
teractions or flows. Each of them will generally produce different orientation 
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Figure 2: Example coverages for two different collector shapes and three different 
a 2 . From left to right a 1 = 1 .0, 0.5 and O.f which according to ^ corresponds to 
order parameter q = O.f 88, 0.668 and 0.982 respectively. The expected direction 
was (p = 7r/2 for round collectors and 7r/4 for square ones. 



probability distributions with variance being typically a nonlinear function of 
experimentally controlled parameters. On the other hand, this study focuses on 
properties of ordered layers regardless of how they were created. Therefore, it is 
worth to introduce separate parameter describing orientational order, irrespec- 
tively of underlying phenomenons promoting parallel alignment of molecules. 
The order parameter can be expressed as follows [||: 



q = 2 



1 



\ " 



N ^ 

i=l 



(xi cos 4> + Di sin < 



(2) 



where n is a number of molecules in a layer, [xi, y.i] describes a unit vector paral- 
lel to orientation of i-th fibrinogen, and <fi denotes a mean direction. Parameter 
q defined as above is for totally disordered layer and rises up to 1 with growing 
orientational order. The maximum value is reached when all molecules inside a 
coverage are parallely aligned. 



3. Results and discussion 

3.1. Coverage ratio 

The maximal random coverage 9 ma x was estimated by measuring the number 
of adsorbed molecules n: 

@max — ^7; • (3) 

Determining its dependence on orientational order was the main aim of this 
paper. Values obtained from simulations were averaged over at least 20 different 
collectors. Results are presented in Figj3] As expected the maximal random 
coverage grows with orientational order, although the overall impact is only 
f0%. The trend of the growth can be successfully fitted with a parabola in 
which the free term is in concordance with the maximal random coverage for 
fibrinogen layers determined in Q. 
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Figure 3: The maximal random coverages dependence on the order parameter 
@. Dots are simulation data whereas the solid line corresponds to square fit: 
Omax = 0.034 • q 2 + 0.295. 



3.2. Autocorrelations and structure 

To study the structure of obtained fibrinogen layers in more detail we anal- 
ysed density autocorrelations and distribution of uncovered places. The density 
autocorrelation function is defined here as: 



G(r) 



2nrp' 



(4) 



where P{r) is a probability of finding two molecules in distance r and p is mean 
density of molecules inside a covering layer. Fig0] shows G(r) function graph 
for three different values of the order parameter. It is noteworthy that in case 




Figure 4: The density autocorrelations for coverages characterised by three 
different order parameters. 



of fibrinogen lack of global order does not mean lack of order at all. It can be 
noticed even on snapshots presented in figH] however, presented G(r) proves 
that. Apart from the statistical noise, all plots look similar. In all the cases 
local maxima and minima are in the same positions and correspond only to 
geometrical properties of the fibrinogen molecule. 
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Another interesting characteristic of a coverage structure is distribution of 
uncovered places. It can be of particular importance when the layer undergoes 
further adsorption of different molecules. To measure this property, the uncov- 
ered places were filled with balls of a different size giving the higher priority to 
larger spheres. It occurs that the size distribution of balls for all values of the 
order parameter q is well approximated by an exponential law: 

N(r) = Aexp(-7r) (5) 

where N(r) is the number of balls having radius r and A is a normalisation 
constant. The decay constant 7 depends on orientational order inside a layer 
as shown in FigO Lower values of 7 for disordered layers mean more large 




Figure 5: The decay constant 7 defined by ((5]) dependence on orientational 
ordering. Dots are simulation data whereas the solid line is a quadratic fit: 
7 = 0.335 • q 2 + 0.049 • q + 0.599. 

uncovered spots, which is in line with the expectations. 
3.3. Kinetics 

There are two factors affecting the adsorption kinetics measured in experi- 
ments. The first one is transport process shifting molecules to a surface proxim- 
ity. It depends on a given experiment's environmental conditions and therefore 
hardly enters the general theoretical analysis. The second factor is varying 
probability of adsorption. It decreases in time due to diminishing area of un- 
covered collector surface. The dependence between adsorption probability and 
temporary coverage ratio is commonly known as Available Surface Function 
(ASF) and it can be easily measured in RSA simulations. Results are shown 
in Fig|51 To compare data for different q, the coverage ratios were normalised: 
6 = 0/9 max . At the limit of low coverage, the ASF(9) is commonly approxi- 
mated by a quadratic fit 0: 

ASF(9) = 1 -C 1 9 + C 2 2 . (6) 

Coefficients C\ and C2 can be calculated using the least squares method. Their 
values for different order parameter q are presented on the inset in Fig|6l Both 
plots in this figure shows that the ASF decays slower for highly ordered layers 
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Figure 6: Available Surface Function for different order parameter q. Inset 
shows coefficients C\ and C 2 used in low coverage approximation ((SJ and their 
dependence on q. Squares and triangles represents data obtained from simula- 
tions whereas lines are exponential fits: C\(q) — 4.709 — 0.015 • exp(4.799g) and 
C 2 {q) = 6.158 - 0.053 • exp(4.412g). 



which is in agreement with the expectations and previous results for maximal 
coverage ratio. 

At the jamming limit, adsorption practically does not depend on transport 
process because its speed is restricted by amount of uncovered spots. Therefore, 
due to it's universality, this case has been extensively studied. For example, for 
spherical molecules it was shown analytically that the coverage's growth obeys 
the so called Feder's law [1, [sjj: 

max - 0(t) ~ r^ d , (7) 

where d is a dimension of a space (collector and adsorbate). Relation ([7]) has 
been proved numerically for a number of dimensions 13, U | and different ad- 



sorbates Q. Since ASF(9) ~ d6/dt, equation ((?]) is equivalent to: 

ASF (6) ~ (1-0)°, (8) 

where a — d + 1. The a values can be obtained from plot |B] usind the least 
squares approximation method. Relation ([5J is valid within 9 — > 9 max limit, 
therefore only data for 9 > 0.8 were used for estimation. Results are shown 
in FigfTJ The statistical errors are quite large due to poor statistics; however 
the correlation coefficient in all cases exceeded 0.9 which means quite good 
agreement with (|5]). As a < 3 the Feder's law is never valid for fibrinogen 
adsorption in 2D. Parameter a fluctuates around 1.5 to rise up to 2.0 for the 
most ordered layer. It suggests that ordered fibrinogen adsorption kinetics is 
similar to one dimensional adsorption. Smaller values of a for disordered phase 
are unexpected and do not agree with comm only used fit for elongated molecules: 
ASF{9) = (1 + ai 9 + a 2 9 2 + a 3 9 3 ){l - 9) A 0. 

4. Summary 

The global orientational ordering inside a fibrinogen layer affects the max- 
imal random coverage ratio, however the overall increase is smaller than 10%. 
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Figure 7: Dependence between a and orientational order q. 

The density autocorrelation does not change as the local orientational ordering 
is present also for globally disordered layers. The analysis of remaining uncov- 
ered spots shows that larger areas are more common in disordered coverages. 
Adsorption kinetics depends on adsorbate alignment and is faster for parallely 
oriented molecules. RSA of fibrinogen does not obey the Feder's law in its orig- 
inal form. For ordered layers the reaction speed is similar to one dimensional 
systems. 

This work was supported by grant MNiSW/0013/H03/2010/70. 
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